CMB Polarization Systematics Due to Beam Asymmetry: 
Impact on Cosmological Birefringence 
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The standard cosmological model is assumed to respect parity symmetry. Under this assumption 
the cross-correlations of the CMB's temperature anisotropy and 'gradient '-like polarization, with 
the 'curl'-like polarization identically vanish over the full sky. However, extensions of the standard 
model which allow for light scalar field or axion coupling to the electromagnetic field, or coupling to 
the Riemann gravitational field-strength, as well as other modifications of field theories, may induce 
a rotation of the CMB polarization plane on cosmological scales and manifest itself as nonvanishing 
TB and EB cross-correlations. Recently, the degree of parity violation (refiected in polarization 
rotation) was constrained using data from BOOMERANG, WMAP and QUAD. Forecasts have 
been made for near-future experiments (e.g. PLANCK) to further constrain parity- and Lorentz- 
violating terms in the fundamental interactions of nature. Here we consider a real-world effect 
induced by a class of telescope beam systematics which can mimic the rotation of polarization plane 
or otherwise induce nonvanishing TB and EB correlations. In particular, adopting the viewpoint 
that the primary target of future experiments will be the inflationary B-mode signal, we assume 
the beam-systematics of the upcoming PLANCK and POLARBEAR experiments are optimized 
towards this goal, and explore the implications of the allowed levels of beam systematics on the 
resulting precision of polarization-rotation measurements. 

PACS numbers: 98.70.Vc 



INTRODUCTION 



Future CMB polarimeters are poised to detect the B- 
mode polarization. It is expected that PLANCK will de- 
tect the lensing-induced B-mode in the near future and 
a host of CMB experiments are expected to detect the 
horizon-scale inflationary B-mode signal if the tensor-to- 
scalar ratio is r « 0.01. This extremely feeble signal is 
susceptible to a variety of foregrounds and systematics 
which may severely degrade the science which CMB ob- 
servations could otherwise do in an ideal- world. Analytic 
expressions have been obtained for several types of beam 
systematics as well as the acceptable levels of beam un- 
certainties, such as differential ellipticity, gain, rotation, 
etc. (e.g. [l]-[4]). We are especially concerned with those 
beam features whose levels and other characteristics are 
not precisely known (e.g., the extent to which the beam 
is elliptical, its gain mismatch, etc) since this cannot be 
remedied in the process of data analysis and will therefore 
add an uncertainty to the CMB power spectra. The level 
of allowed beam uncertainty is determined mainly by the 
type of science one could do with B-mode polarization; 
these are associated with inflation and lensing-related sci- 
ence such as neutrino masses, dark energy equation of 
state, etc. It might be possible that an actual beam cal- 
ibration will result in a better calibration than the mini- 
mum requirements for the science. As it turns out, for as 
low a tensor-to-scalar ratio as r = 0.01 the most stringent 
experimental requirements on the allowed level of beam 



uncertainty come from the inflationary science. Not only 
will B-mode observations open a unique window to infla- 
tion and gravitational lensing, it also offers a rare oppor- 
tunity to peer into two of the most fundamental under- 
lying symmetries of physics, namely - parity and Lorentz 
invariance. With an horizon-scale optical path we have 
the longest possible leverage-arm to detect the accumu- 
lated effect of 'cosmic birefringence' (CB), a continuous 
rotation of the polarization plane of the electromagnetic 
radiation. Lorentz- and/or parity- violating terms in the 
electromagnetic or gravitational sectors of a hypothetical 
unified model of the fundamental interactions of nature 
may induce such a CB (e.g. [5]- [10]). Searching evidence 
for CB in the CMB power spectra have been the objec- 
tive of several recent works ([11]-[16]) which analyzed the 
data in an ideal, systematics-free, world. 

The purpose of this work is to highlight the importance 
of controlling beam systematics, mainly their effect on 
the polarization state of the CMB and in the context of 
CB. A tantalizing idea is to use the TB and EB cross- 
correlations in the CMB to detect Lorentz-symmetry 
breaking or parity-odd features. However, these cross- 
correlations are especially susceptible to even the small- 
est degradation in the B-mode (due to temperature leak- 
age or beam rotation) and therefore a careful study of 
the implications beam systematics may have on the CB 
science is well motivated. 

The outline of this paper is as follows. Section 2 sum- 
marizes CMB power spectra and relevant Fisher matrix 
relations. The effect of beam systematics is discussed in 
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section 3. We describe our analysis and results in section 
4 followed by a discussion in section 5. We conclude in 
section 6. Our basic calculations are supplemented with 
two appendices: In Appendix A we estimate the resulting 
level of beam ellipticity assuming a standard procedure 
of telescope beam calibration (which determines the level 
of TB and EB correlations induced by beam ellipticity) . 
Appendix B is dedicated to a brief discussion of the ran- 
domness of beam-pair orientation in the focal plane and 
its implications on the spurious power spectra. 

CMB POWER SPECTRA AND CB 

Before we describe the effect of non-standard physics 
on the CMB, let us first review the basics of the CMB 
statistics in the standard cosmological model because de- 
tecting this physics relies on deviations from the standard 
statistics of the CMB. The temperature T and the Stokes 
parameters Q &z U are expanded in scalar and rank ±2 
tensor spherical harmonics 

T(n) = ^aLi;m(n) 
{Q±iU){h) = '^a±2,im ±2Yim{n) (1) 

where the E and B (scalar and pseudo-scalar, respec- 
tively) modes of polarization have expansion coefficients 
Elm and Bim which can be expressed in terms of the spin 
±2 coefficients 

a±2,lm = afm^ O'hn- (2) 

It is assumed, and well motivated by cosmological infla- 
tion, that the temperature anisotropy and polarization 
modes can be drawn from gaussian distributions fully 
characterized by vanishing means and standard devia- 
tions described by the power spectra. This implies that 
the ensemble averages of the expansion coefficients vanish 

{afj = (3) 

where X can assume cither T, E or B and that the widths 
of the gaussians satisfy (by the assumption of statistical 
isotropy which is the multipole-space dual of real space 
isotropy) 

{O'Vm'O'hn) = Cf^^ Sll'Smm' (4) 

and 

m=—l 

In the standard model E and B are pure parity states 
(even and odd, respectively). The correlation, over the 



full sky, of the B-modc with either the temperature or E- 
mode polarization vanishes in this case. However, if the 
polarization plane is rotated for some reason (cosmolog- 
ical or otherwise), there will be a E-B mixing which will 
induce the 'forbidden' TB and EB power spectra (leakage 
from the underlying TE and EE cross-correlations, re- 
spectively). Let us see how this effect takes place. While 
rotation of polarization leaves the temperature invariant, 
the polarization state will be affected. Consider for exam- 
ple the flat- sky approximation, which holds for small an- 
gular scales only (and is fully applicable on beam scales), 
but enables a transparent derivation of the following rela- 
tions between the original and observed degraded power 
spectra. In this approximation the Stokes Q and U pa- 
rameters, which fully describe the polarization state of a 
linearly polarized radiation, can be expanded in Fourier 
space as follows 

{Q ± if/)(x) = ^ f{Ei± iBi)e*' "e±2i(*'-'^^)(i2i_ ^g) 
27r J 

Q ±iU are pure spin ±2 states and under a rotation by 
—7 pick up phases e^'^^'^ . As seen from Eq.(6) this results 
in an E-B mixing 

E' = cos 27 -B sin 27 

B' = £; sin 27 -I- B cos 27 (7) 

and correlating with T' = T we obtain the following 



power spectra 
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Note that no assumption has been made here as to the 
origin of this rotation, namely whether or not it is of cos- 
mological origin. In the literature 7 is identified with the 
CB rotation angle (see e.g. [15] and [16]). We will see in 
the next section that, in practice, one has to account for 
other, non-cosmological effects which mimic cosmologi- 
cal polarization-rotation and may bias the results for the 
best-fit 7 if not properly accounted for (or, to a lesser 
extent, change the uncertainty of the inferred rotation 
itself). The main feature of Eq.(8) is that the 'forbid- 
den' C^^ and C^^ do not vanish in general, and it is 
this property of the power spectra which will be used (in 
the near-future, once we have high signal-to-noise mea- 
surements of TB and EB correlations) to infer the level 
of CB. Current weak upper limits on 7 are derived from 
comparing Cf^ with Cf^ and Cf'^ , all of which are pri- 
marily sourced by density perturbations. In the absence 
of CB, Cf^ directly fixes Cf^ and Cf ^ and the latter 
add little new cosmological information in the standard 
model with no reionization. Any deficit in the latter is 
ascribed to E-B rotation. 
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THE EFFECT OF BEAM SYSTEMATICS 

No CMB observation is perfect and is always con- 
fronted, at some level, by optical beam systematics. 
Beam mismatch can induce temperature leakage to po- 
larization as well as mode-mixing between the E- and 
B-modes ([1], [17]). Temperature leakage to polariza- 
tion is considered the most pernicious systematic due to 
the fact that the primordial temperature anisotropy is at 
least an order of magnitude larger than the E-mode and 
two-three orders of magnitude larger than the B-mode. 
Therefore, even a small level of temperature leakage to 
polarization significantly contaminates the signal. For 
CB science, however, we are most concerned with cross- 
polarization and we open this discussion with the effect 
of beam rotation which mixes E and B-modes ([1], [17]). 
Beam rotation by an angle e trivially biases 7 at this level 
which implies that the uncertainty on the inferred 7 can- 
not be smaller than the bias induced by the uncertainty 
in beam orientation, e. Threshold values for e have been 
found for POLARBEAR, PLANCK and CMBPOL [3] 
based on the requirement that both inflation and lensing 
science are not degraded above certain level. However, 
this requirement does not necessarily protect the CB sci- 
ence from beam systematics. 

CMB polarimetcrs such as PLANCK will have suffi- 
ciently high instrumental sensitivity to allow a statisti- 
cally significant detection of C^^ (in a systematics-free 
world), the 'smoking gim' for CB. As will be shown be- 
low, PLANCK and POLARBEAR have the raw sensitiv- 
ity to detect rotation angles as small as few arcminutes. 
Any beam rotation by a larger angle will wash out the 
cosmological rotation of the polarization plane if CB is 
small enough. In these high-sensitivity experiments high 
polarimetric fidelity is required. 

In the following we briefly discuss what impact two 
other beam systematics have on CB. These two effects, 
due to beam ellipticity and pointing error, leak temper- 
ature to polarization and are not a simple mixing of 
polarization as is the case with beam rotation. Since 
these two offsets will have a unique spectral-dependence 
in Fourier space (the systematic B-mode power spectra 
due to differential pointing and ellipticity scale as oc l^Cf 
and oc l^Cf , respectively) they can be statistically dis- 
tinguished and separated (with some residual 'noise'), in 
principle, from CB which inherits the spectral shape of 
Cf^ Eqs.(7) will not be preserved in this case and we do 
expect that the likelihood function will change its shape 
and not be only biased. Table I summarizes the effect 
of differential rotation, pointing and ellipticity studied 
in [17] on Cf^ and Cf'^ . It should be mentioned here 
that two other beam systematics considered in [17] and 
[3], differential gain and beamwidth, do not contribute 
to Cf^ and Cf ^. They do, however, affect Cf^, Cf ^ 
and Cj^^ , but since our main goal here is to test the 



susceptibility of Cf^ and Cf"^ to beam systematics and 
the implications it has on the inferred 7 we do not ad- 
dress the effects of differential gain and beamwidth in 
this work. 



ANALYSIS AND RESULTS 

An easy way to assess parameter uncertainty (7, the 
CB rotation angle in our case) is to invoke the Fisher 
information-matrix approach (though not in the presence 
of signiflcant bias, as discussed below). The elements of 
the Fisher matrix are defined as follows 

1 pir' 
F^J = 2 + l)/afe.'n:ace[C-^ — ] (9) 

where the 4-D symmetric covariance matrix C; (per each 
multipole /) is defined as 
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cr 




cr 
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Cfd J 



(10) 



where Cf''' and Cr refer to the power spectra associ- 
ated with the lensing deflection angle [3]. The partial 
derivatives are taken with respect to the parameters Aj. 
Here we allow nonvanishing TB and EB correlations and 
include lensing extraction. The dimensionality of the 
Fisher matrix is N x N where N is the dimension of 
parameter space of the cosmological model in question. 
The la error on the parameter A, is given by 



(11) 



The la error obtained by this procedure only sets a lower 
limit on the actual error. In general, when parameter 
degeneracy is high, as well as in the presence of significant 
bias, one resorts to Monte Carlo Markov Chain (MCMC) 
analysis [3]. The terms Cr ^nd Cr drop from the 
matrix C (Eq. 10) in the standard model, but in general, 
especially here, they do not. 

A small bias in a parameter Xi compared to the statis- 
tical uncertainty, is given by [2,3] 



AA, = (Af«)-(Ar-) = ^(F-i)..B,. 



where the bias vector B can be written as 



f)f^cmb 



(12) 



(13) 



^ij = cov(Q,Cf) and C; is a vector containing all six 

power spectra. Our analysis includes lensing reconstruc- 
tion and is essentially a generalization of [3] . 
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As pointed above already, the underlying assumption 
of this work is that the primary target of upcoming 
CMB experiments is the inflationary B-mode; lensing- 
induced B-modc as well as the B-mode polarization in- 
duced by the Chern-Simons-type interaction term are 
only secondary. In addition, it was shown already [3] 
that the inflationary signal is more prone to beam sys- 
tematics than the lensing signal is and therefore it was 
assumed that optimizing the beam systematics to the 
former should be sufficient for the latter. Here, we use 
the analytic expressions for beam systematics found in 
[17] to determine the resulting bias in the inferred CB 
for both PLANCK and POLARBEAR. Before present- 
ing the results of our Fisher-matrix analysis it is con- 
structive to plot the systematic Cf^ and Cf^ . These 
are shown for both PLANCK and POLARBEAR in Fig- 
ure 1. Black (dotted), blue (dashed) and yellow (dotted- 
dashed) curves correspond to the systematics induced by 
differential pointing, ellipticity and beam-rotation, re- 
spectively (note that the effect of pointing on the EB 
cross-correlation is negligible and falls below the scale 
shown). For all the effects we assumed worst-case scan- 
ning strategy and polarimeter orientation as will be dis- 
cussed below and extensively discussed in [17] and [3]. 
Red (continuous) curves show levels of Cf^ and Cf ^ 
which these experiments could have detected (obtained 
from the Fisher-matrix analysis) in the absence of sys- 
tematics. We refer to them as 'nominal CB' detec- 
tion. For comparison, the green (continuous) curve corre- 
sponds to a much larger rotation of the polarization plane 
that will result in Cf^ and Cf^ roughly three times 
larger than the largest systcmatics-induced TB and EB 
correlations. It is evident from the figure that realizing 
the nominal potential of PLANCK and POLARBEAR to 
detect CB via the TB and EB correlations will require a 
much better control of systematics than the level set by 
the inflationary-induced B-mode detection requirements, 
for example. 

The results of our Fisher-matrix analysis are shown in 
Table II. Nominal uncertainties in 7 (i.e. the statisti- 
cal error in a systematics-free experiment) are compared 
with the bias induced by differential pointing and ellip- 
ticity, A7. The values for the differential pointing p and 
ellipticity e are adopted from [3], consistent with the as- 
sumption that PLANCK and POLARBEAR are opti- 
mized for inflation and lensing science. In the analysis 
we used all information from multipoles I < 1200. Figure 
1 suggests that the bias decreases if Imax is smaller but 
this of course comes on the expense of statistical uncer- 
tainty. We repeated the analysis with lower Imax which 
indeed resulted in lower bias but even for Imax 

as small 

as 200 it was Lmacceptably large. The effect of beam 
rotation results in a pure bias in the inferred 7 and is 
therefore trivial and not included in the table. Again, 
we find that beam calibration levels that suffice for the 
inflationary and lensing science are unfortunately insuf- 



flcient for the CB science and this is the main conclusion 
of this paper. Especially worrisome in this context are 
beam rotation and differential ellipticity. Both system- 
atics arc independent of the scanning strategy and there- 
fore cannot be easily mitigated. In addition, beam ro- 
tation induces both EB and TB correlations which have 
the same spectral shape in multipolc-spacc as the CB; 
one cannot distinguish a beam rotation from CB. While 
it is clear from our analysis that the allowed ellipticity 
based on the inflation+lensing requirements may not be 
sufficient for CB purposes, beam calibration in a real ex- 
periment will result in a better control of ellipticity than 
the minimum allowed by inflation/lensing. We show in 
Appendix A that both PLANCK and POLARBEAR wiU 
benefit from a reasonable control of ellipticity which will 
allow beam ellipticities lower than the minimum required 
for inflationary and lensing science and will therefore re- 
sult in an essentially unbiased CB detection from the 
E-B correlations of the CMB in case the nominal val- 
ues of the respective 7 are considered. However, TB is 
still prone, even in this ideal case, to beam ellipticity. 
Beam rotation remains the paramount concern, and of 
all types of systematics we identify it as the main ob- 
stacle towards CB detection; the nominal 7 values are 
typically a factor « 10 smaller than the allowed beam 
rotation (Table 2). Recently, it was shown how optimal 
estimators, which can filter spatially-dependent rotation, 
can be used to filter-out a spatially-dependent rotation 
due to the non-standard statistics it induces by inducing 
I — I mode coupling [18]. However, the constant beam 
rotation we consider here cannot be distinguished from 
CB by this method. 

The systematics induced by differential pointing de- 
pend, to first order, on the scanning strategy and for a 
non-ideal, yet uniform, scanning strategy we encapsulate 
the relevant information in the two scanning functions /2 
and /s ([17]) whose exact definitions are irrelevant to this 
discussion since we always assume the worst case scenario 
([3])/i = 27r. 

The cosmological model adopted in our numerical cal- 
culations closely follows the WMAP5 results [15]; the 
baryon, cold dark mater, and neutrino physical energy 
densities in critical density units Q.i,h? = 0.021, flch^ = 
0.111, riiyh"^ = 0.006. The latter is equivalent to a total 
neutrino mass M^, = Yl'f^i m^^i =0.56eV. Dark energy 
makes up the rest of the energy required for a spatially- 
flat universe. The Hubble constant, dark energy equation 
of state and helium fraction are, respectively, Hq = 70 
km sec~^Mpc~^, w = —1 and Yhb = 0.24. h is the 
Hubble constant in 100 km/scc/Mpc units. The optical 
depth to reionization and its redshift are r^e = 0.073 and 
Zre = 12. The normalization of the primordial power 
spectrum was set to As = 2.4 x 10~® and its power law 
index is Ug = 0.947. 
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DISCUSSION 



CMB observations, especially those which will have the 
sensitivity and fidelity to detect the 'forbidden' power 
spectra C^^ and C^^ , will be able to set tight con- 
straints on CB, a phenomena related to the bedrock of 
fundamental physics and its symmetries. Cosmology im- 
proves on terrestrial experiments in this context due to 
the long optical path lengths which enable the detec- 
tion, in principle, of extremely small effects. However, 
since a statistically significant detection of CB heavily 
relics on measuring the TB and EB 'smoking gun' corre- 
lations, and due to the fact that these involve the suh-^K 
B-mode signal which is prone to numerous systematics, 
a credible detection of non-vanishing 7 should account 
for these possible sources of confusion. We saw that 
beam-rotation, differential ellipticity, as well as differ- 
ential pointing, gain and beamwidth (to a lesser extent) 
can bias the inferred 7 but also change the uncertainty 
(in the case of differential pointing and ellipticity) due to 
their direct effect on the TB and EB correlations. While 
the differential pointing is very small due to its shal- 
low 1-dependence (compared to other systematics con- 
sidered here) and the effect of beam ellipticity can be 
partially harnessed (Appendix A), it is the beam rota- 
tion effect which mainly contaminates CB. In particular, 
we have shown that CMB experiments which are opti- 
mized for inflationary and lensing science may not be 
adequate for the CB science that requires a much bet- 
ter control of beam rotation. The nominal 7-detcction 
with PLANCK and POLARBEAR is 3.8 and 1.1 ar- 
cminutes, respectively. Beam rotation should therefore 
be controlled to 0.75 and 0.21 arcmimites if a 5(t detec- 
tion of CB is required. This systematic effect cannot be 
mitigated by scanning strategy. The effect of differential 
pointing (which is subdominant to the effect of differ- 
ential ellipticity) can be further suppressed by scanning 
strategy mitigation [17]. The differential-ellipticity, to- 
gether with the beam rotation effect, seem the most per- 
nicious for TB spectra (Fig. 1). We consistently assumed 
the worst-case-scenario as we did when we employed the 
constraints on C^^ by assuming the tilt of the polariza- 
tion direction to the ellipse major axis is 45'' [3] (see also 
Figure 2). However, an estimate of the real uncertainty 
on beam ellipticity (Appendix A) shows that, in prac- 
tice, beam ellipticity should not be a major problem to 
TB measurements by PLANCK and POLARBEAR. The 
EB cross-correlations, however, should vanish if we take 
ip = 45*^ = 9. Nevertheless, for the purpose of illustra- 
tion, and in the plots only, we assumed ip = 22.5° = ^ so 
as to maximize the EB correlations for the given beam 
parameters [3]. We further explored the effect of ran- 
domness of the angle ^ in a multi-pixel elliptical-beam 
experiment (Appendix B) and found that our worst-casc- 
scenario is actually representative of an average multi- 



pixel experiment. Similar considerations were also ap- 
plied to beam rotation. 



CONCLUSION 

CB is an interesting ancillary science that future CMB 

polarimcters will target. We considered in detail two such 
experiments: PLANCK and POLARBEAR. Our past ex- 
perience taught us that controlling beam systematics to 
a sufRciently high precision as to aim at the B-modc de- 
tection (in case the tensor-to-scalar ratio is r ~ 0.01) 
will also automatically gurantee a high fidelity lensing- 
induced B-modc measurement. However, as we show 
here, this minimum requirement will in general not suffice 
for a credible CB detection via the 'forbidden' TB and 
EB correlations. The likely reason for this is that while 
BB correlations are quadratic in the beam-imperfection 
parameters (e.g. oc e^,e^) the TB correlations are only 
linear in these small beam parameters (and the EB cor- 
relation is in general noisy) . This implies that for a given 
ellipticity or beam rotation the fractional bias in the TB 
cross-correlations will be larger than the the correspond- 
ing fractional bias in the BB power spectrum by 0(l/e) 
and 0(l/e), respectively. 

Whereas the effect of differential pointing on the TB 
and EB power spectra, we showed, is negligible, the ef- 
fect of ellipticity is larger and can compromise the CB 
science if only the minimal requirement of the infla- 
tion/lensing science is satisfied. However, we demon- 
strated (Appendix A) that in practice the expected level 
of ellipticity uncertainty is relatively small (if the beam 
is standardly calibrated against a bright point source). 
The most pernicious effect is due to beam rotation which 
precisely mimics BC, has the same ^-dependence, and 
cannot be mitigated by idealizing the scanning strategy. 
Our conclusion is that while TB and EB are unique indi- 
cators for new-physics in principle, they can in practice 
be excited by imperfect beams and in order to realize 
the promising potential of the high-sensitivity and fine- 
resolution PLANCK and POLARBEAR at CB detection 
beam rotation has to be controlled to the sub-arcminute 
level. Unless beam rotation is controlled to the arcminute 
level, a conservative approach, which docs not use Cf^ 
and C^^ , may be more adequate for the CB science. 
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FIG. 1: TB and EB systematics compared to the nominal CB signals (i.e. those we expect to obtain in the absence of beam 
systematics) : black (dotted), blue (dashed) and yellow (dotted-dashed) curves refer to the systematics induced by differential 
pointing and ellipticity and by beam rotation, respectively. Red and green (continuous) curves represent the cross-correlations 
in the nominal and realistic CB cases, respectively. 'Realistic' is defined as the CB that would be roughly three times larger 
than the largest systematic. Dashed-dotted curves represent the absolute values of negative T-B cross-correlations in certain 
/-ranges. All signals are calculated based on the assumption that the beam systematics of PLANCK and POLARBEAR are 
optimized for the detection of the inflationary signal at the r — 0.01 level [3]. These are used in the analytic expressions of 
[17] for Cf^ and Cf"^. It is clear from the upper panel that the main contaminant is due to differential ellipticity and beam 
rotation. 



Appendix A: Forecasted Beam Ellipticity 

It is instructive to estimate the expected uncertainty of 
beam ellipticity which may resuh from a standard beam- 
calibration procedure with a nearly black-body point- 
source such as Jupiter (Tp « 150^, 9p ~ 0.5 arcmin). 
By Wiener filtering [4] a map of the expected signal one 
expects to recover the source image with a signal-to-noise 
level 

\n) -J P(l) (2^)2- 



(A.l) 

where S{1) is the Fourier transform of the point source 
and P{1) is the instrumental noise, i.e. the (S/N)'^ is 
the ratio of the signal and noise power-spectra integrated 
over all accessible multipoles in the experiment in ques- 
tion. The higher S /N the smaller are the uncertainties in 
the recovered beam parameters. The Fourier transform 
of the convolved calibration source is 

(A.2) 
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where we assume an elliptical gaussian beam with prin- 
cipal axes (Tx and ay, and pointing p. This results in 

fsY_ 2in(2) fTA'feA\,. 

(A.3) 

where rj is an experiment-specific optical-efficiency pa- 
rameter (< 1) and ft is the fraction of observation time 
dedicated to beam calibration (we adopt the value f /30) . 
Afc is the instrumental temperature equivalent noise, 9h 
is the beamsize and e is its ellipticity. Since the point- 
ing merely adds a phase to the beam function it drops 
from the expression for S/N. Similarly, S/N is also in- 
dependent of the beam rotation angle since temperature 
measurements are insensitive to e. Calibrating the beam 
rotation with a polarized source will not work here since 
by analogy to Eq. (A.l) the rotation only adds a phase. 
Therefore, the following procedure, which is based on 
S/N considerations, will be used to determine the un- 
certainty of e only. To estimate these uncertainties we 
require that varying the beam parameters results in sig- 
nal changes smaller than the noise. That is, we allow the 
uncertainty in the beam's ellipticity to grow from until 
the (S/N)"^ changes by unity 

{S/Nf ^ {S/Nf + 1. 

(A.4) 

This condition readily yields the uncertainty in beam pa- 
rameters 

(A.5) 

Assuming Jupiter is a 150^ black body with a char- 
acteristic scale of Op — 0.5' , typical beamsizes 6t are 
5' (4') and instrumental noise Af, l.TjiK (lO/iX) for 
PLANCK (POLARBEAR) as well as 5% optical effi- 
ciency rj we obtain the following estimates on expected 
uncertainty in beam elhpticity e = 1.1 x 10"'^ (PLANCK) 
and e = 1.2 X 10"'' (POLARBEAR). Using the scahng 
of ellipticity-induced TB and EB correlations (oc e and 
oc e^, respectively), it is clear that based on this estimate 
beam ellipticity is unlikely to be the dominant contami- 
nant of the EB correlations but the TB correlations will 
still be somewhat biased even for the optimal elliptic- 
ity we find here (0(10-^) for PLANCK and 0(10^"*) for 
POLARBEAR). 

Appendix B: Focal Plane Considerations 

Previous works [1, 2, 3, 17] assumed, for simplicity, 
single-pixel experiments. However, in practice the beam- 
pairs are scattered in the focal plane with different polar- 
ization orientations and other beam parameters such as 



gain, beamwidth, ellipticity, etc. We briefly discuss how 
the random orientation of beam polarization only mildly 
suppress spurious B-mode polarization. A schematic fo- 
cal plane is depicted in Fig. 2. We assume all beams 
have the same ellipticity for the purpose of illustration 
but they differ by their orientations with respect to the 
radius- vector which connects the center of the focal plane 
to the beam center (denoted x) and the polarization axis 
makes an angle if) with the ellipse major axis. We further 
assume that the major axes of the beams are nearly or- 
thogonal (xi — X2 ~ 90°) and similarly their respective 
polarization directions satisfy (V'l — ~ 90°), where the 
subscripts stand for 'beam 1' and 'beam 2' of each pair, 
respectively. We calculate the covariance of the individ- 
ual systematic power spectra 

cov{cr,cj)^6iA{{crf)~{{cr)f] 

(B.l) 

where Y stands for cither BB, TB or EB and the an- 
gular brackets denote averages over the angles ^ and x- 
Since the number of pairs is finite, all the estimates we 
obtain here for the suppression of systematics by assum- 
ing infinitely many detector pairs constitute only upper 
limits. Nevertheless, we show that even if the number of 
detectors is assumed infinite, the gain in systematic sup- 
pression is insignificant. Taking the ratio of the covari- 
ance to the worst-case Cf , Cf^ and Cf ^ (^/; -f ;x = 45° , 
45° and 22.5°, respectively) it is straight forward to show 
that with an infinite number of pairs the systematic Cj^ , 
Cf^ and Cf ^ drop to 71%, 71% and 50% of their worst- 
case values - an insignificantly small change to our con- 
clusions. We can apply a similar calculation to beam 
rotation, in this case averaging is carried out over the 
angle e (random orientations of the pair subject to the 
constraint that the polarization directions of the two po- 
larimeters remain orthogonal). Here we obtain that Cj^ , 
Cf^ and Cf ^ are suppressed to 35%, 71% and 50% of 
their worst-case values. Again, this will not quanlita- 
tively change our conclusions. 
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effect 


parameter 


Acr 




gain 


a 








monopole 


M 








pointing 


P 


^[1 + jQilp)]J2ilp)seCf 
+^seJl(lp)CT 


SgCgJi{lp)Cf 


quadrupole 


e 


-Io{z)Ii{z)s^CT 




rotation 


e 


sm{2e)Cf'^ 


isin(4e)Cf 



TABLE I: The contribution of the systematic effects to the 
power spectra Cf^^ , Cp^ assuming the underlying sky is not 
polarized (except for the rotation signal when we assume E-, 
and B-mode polarization are present) and ideal sky scanning. 
The definitions of sg and as well as /a can be found in [17] 
where similar expressions for Cf^, Cf^ and C^^ are also 
found and which we employed in this work. 



Experiment 


^nominal [jgg] 


parameter 


A7[deg] 


POLARBEAR 


0.0179 


P 


0.0094 


e 


56.07 


PLANCK 


0.0641 


P 


0.0060 


e 


93.9 



TABLE II: Estimated bias in the inferred rotation angle as- 
suming POLARBEAR and PLANCK were optimized to de- 
tect the inflationary B-mode signal associated with r=0.01 
([3]). The bias induced by differential ellipticity is much larger 
than the naive uncertainty, highlighting the need in careful 
data analysis when interpretation of polarization plane ro- 
tation is required. The relatively small bias induced by the 
pointing is negligible compared to the ellipticity-induced bias 
and this is consistent with the findings of [3] and explained 
by the steep 1-dependence of the ellipticity systematics. 



FIG. 2: Layout of typical beam pairs in the focal plane: dashed lines represent radius- vectors from the center of the focal plane 
through the center of each pair. Continuous lines represent the ellipse major axes and the angle the major axis makes with the 
radius-vector is x- The polarization axis (dotted line) makes an angle tp with the major axis. In [17] and [3] it was assumed 
that X = 0: therefore every in our analytic expressions (Table 1) should be replaced with ip + x- 



